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Abstract. The Discrete Absorption Components (DACs) commonly observed in the ultraviolet lines of hot stars 
have previously been modelled by dynamical simulations of Corotating Interaction Regions (CIRs) in their line- 
driven stellar winds. Here we apply the kinematic DAC inversion method of Brown et al. to the hydrodynamical 
CIR models and test the reliability of the results obtained. We conclude that the inversion method is able to recover 
valuable information on the velocity structure of the mean wind and to trace movement of velocity plateaux in the 
hydrodynamical data, though the recovered density profile of the stream is correct only very near to the stellar 
surface. 
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1. Introduction 

The time-averaged structure of line-driven winds of 
hot OB stars can be adequately described by time- 
independent wind models (eg. Hillier & Miller [TI)98I Vink 
et al. I^nnni Pauldrach et al. PUTO . However, winds of 
hot stars are well-known for the variability manifested in 
different domains of the electromagnetic spectrum. One 
of the most striking features in the spectra of hot star 
winds is the presence of Discrete Absorption Components 
(DACs). These DACs are superimposed on the (unsatu- 
rated) absorption part of P Cygni profiles and move from 
the lower blue-shift to higher blue-shift part of P Cygni 
profiles (e.g. Prinja & Howarth H^M Massa et al. [f^ 
Kaper et al. lTMTl Prinja et al. lMT^ . 

The main DAC features, i.e. their recurrence and slow 
acceleration, can be explained by the radiatively driven 
hydrodynamical simulations of Cranmer and Owocki 
(|1996|l . They assumed that DACs arise from Corotating 
Interaction Regions (CIRs), as first suggested by MuUan 
P984(l . The origin of CIR formation is not yet known, 
though wind modulations by stellar pulsations or mag- 
netic fields are often cited as possible drivers (Kaper et al. 
[TW71 de Jong et al. 120011) • Thus it is worthwhile to study 
the DAC phenomenon from data diagnostically, indepen- 
dent of any underlying driver model. A first step in this 
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direction was made by Brown et al. H2(JU3fl . who devel- 
oped a kinematic data inversion method for the inference 
of absorbing stream properties from data on recurrent 
DACs. The optical depth variations calculated using the 
hydrodynamical structures of the Cranmer and Owocki 
l|1996p models on the other hand arise mainly due to the 
azimuthally localized velocity plateau (so-called Abbott 
kinks) propagating upstream in the stellar wind. These 
velocity plateaux are caused by the presence of radiative- 
acoustic waves (the so-called Abbott waves, Abbott ri980|l 
moving with velocity higher than the sound speed. Thus, 
it is interesting to study how well the kinematic inver- 
sion method for line-profiles performs using the profiles 
from hydrodynamical CIR model calculations, specifically 
to see what the method recovers in terms of density and 
velocity variations. 



2. Kinematic inversion formulation 

The basic method of the DAC data inversion of Brown 
et al. H2003|l can be summarised as follows. They assume 
that DACs arise as the result of an azimuthally varying 
mass loss rate at the stellar surface, giving a density pat- 
tern in the wind that also varies with azimuth, and is then 
transported through the wind along the streamlines of an 
axisymmetric velocity law Vr = Vr{r). The assumed ax- 
isymmetry of the wind velocity has several consequences: 
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1. From mass continuity it follows that the radial mass 
flux nVr falls off with radius along any streamline 
as l/r^, which is the purely geometrical factor ex- 
pected for essentially spherical wind expansion (which 
we view only in the equatorial plane - see below) . 

2. From^ it further follows that the density itself at any 
point along a streamline is related to the mass flux uoVq 
at the point on the surface from which that streamline 
emanates by 



n{r) 



voR^ 

2""-0 



(along a streamline), 



(1) 



where R is the stellar radius. 
3. The streamlines from any two points on the surface 
(at 01 and 02, say) will always have identical shapes 
for an axisymmctric velocity law, but will be shifted 
relatively in azimuth by 4>2 — 4'i ^-t all radii. This means 
that the density pattern introduced into the wind at 
the surface will be preserved at all radii, but will be 
shifted in </> by the change in azimuth A0(r) of the 
streamlines due to rotation. In the corotating frame 
this is given by 



A0(r) 



(2) 



where VL{r)r — is the rotational speed in the static 
frame and flo — fl{R) is the angular velocity of the 
surface. (The density will also be scaled according to 
Eq. iQJ along the streamlines.) 

Putting these points together we can relate the physical 
density n(r, (j)) of the wind at any point to the surface 
density profile no{if)) = n{R, ip) by 



(3) 



so that the optical depth variation at the velocity v ~ 
Vr{r) corresponding to this radius is (from r = n/v'^, in 
suitable units) 



easy to determine P(w) and A0(r). Integration of Eq. (0J 
over (j) at fixed Vr gives 



2tt 



V{v) oc / t{v, (f) I 



(6) 



Alternatively, if it is possible to identify the path of a 
single matter stream in the {vr, 0)-plane, such as the peak 
of the density pattern in azimuth at each t;^., the optical 
depth variation along this matter stream can be used to 
define V, although this is likely to be more sensitive to 
data errors on the dynamical spectrum. To obtain A0(r) 
it is necessary to identify just such a matter stream, and 
to define its path in the [vr, (f)) plane. 

3. Dynamical DAC prediction and kinematic 
model inversions 

3.1. DAC prediction 

For the hydrodynamical prediction of DACs we apply a 
similar method to Cranmer & Owocki H1996|l . We assume 
that DAC features in the spectra of OB stars are caused 
by CIRs in the stellar wind induced by a hot spot on the 
stellar surface. Although CIR formation in a radiatively 
driven wind is inherently three-dimensional, to make the 
problem more tractable we confine the hydrodynamical 
calculations only to the equatorial plane, — 7r/2. 

In order to obtain static DAC models we relax the 
time-dependent hydrodynamical equations to a steady- 
state. For the solution of time-dependent hydrodynami- 
cal equations wc use the code VH-l developed by J. M. 
Blondin and collaborators. The equations solved are equa- 
tion of continuity (conservation of mass) in the form 



dp 1 / 



1 



d 



at 



(pv^) = 0, 



(7) 



and the conservation of radial (r) and azimuthal ((/)) com- 
ponents of momentum, 



T(i;„0)^n(r,0)/V = 7'(i;,.)"o(0-A0(r)), (4) dv^ ^ ^^dv^ ^ ^^dv^ ^ ^J> 



1 dP 



where 



vqR^ 



(5) 



gives the overall (^-independent) weighting of the optical 
depth at each velocity and v'^ is the velocity gradient. 

The inversion method of Brown et al. (2003) uses the 
variation of the observed optical depth (using the Sobolev 
approximation Eq.Q) with time (i.e., with observer az- 
imuth, (f> = flot) to determine the optical depth weighting 
function V{vr) and the spiral law A0(r), which then give 
the wind velocity law Vr{r) and rotation law r2(r). Finally, 
from the known optical depth variations with azimuth at 
fixed radius (i.e., velocity) it is possible to calculate the 
surface density profile no((/)o). 

When the observed dynamical spectrum is consistent 
with the assumption of an axisymmetric velocity law it is 



dt 



7r+9os + 

r p or 



^ lines 



dr rsm 

dv^ f0 dv^ _ VrV4, _ 1 dP 
dr r sin 9 dcj) r pr sin 9 d<j) ' 



(8) 

(9) 



where p is the local wind density, Vr and are the radial 
and azimuthal components of wind velocity, t is time, and 
gas pressure P is evaluated using the perfect gas equation 
of state. 

The effective gravitational acceleration is given by 



ffoff 



GM{1 - F) 



(10) 



where G is the gravitational constant, M is the stellar 
mass and the Eddington factor F accounts for the radia- 
tive force due to electron scattering. The radiative accel- 
eration due to fines in the Castor, Abbott & Klein H1975|l 
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approximation is given by 



„lincs 



GMT ( n^/W 



IQi 



dVr 



dr 



(11) 



where fc, a and 5 arc line-force parameters f Abbott 119821 
Puis et al. 12000,) . / is the finite-disk correction factor 
(Friend & Abbott EHHl Pauldrach et aL lB??^ . is the 
electron number density, W is the dilution factor, Ko is 
the opacity due to the electron scattering, and wth is a 
fiducial ion thermal speed (by definition hydrogenic) . The 
function A{i]j) accounts for the force enhancement due to 
the spot, where the corotating azimuthal angle is given 

by 

^ = (j)~n^t, (12) 

with the star's rotational angular velocity fig = V-^ot/R*- 
The force amplitude has the form 

A(^) = 1 + Aoexp [-(^-i^o)V^'] , (13) 

where Aq is the dimensionless spot amplitude (correspond- 
ing to the flux enhancement), ipo is the corotating az- 
imuthal position of the center of the spot and a describes 
the spatial extent of the spot on the star. Thus, contrary 
to the paper of Cranmer & Owocki 1)1996(1 . we have ne- 
glected limb-darkening. 

The flow variables are specifled on a fixed two- 
dimensional mesh in radius and azimuth. We use 250 ra- 
dial zones, extending from r = i?^ to r = Hi?* with higher 
concentration of zones near the stellar base. The spatial 
zone width radially increases by 1% per zone. The az- 
imuthal mesh contains 180 constantly spaced zones, rang- 
ing from to 27r. The boundary conditions are specified in 
the same way as in Owocki et al. p994|l and Cranmer & 
Owocki LiggSJ . As an initial condition we use the relaxed 
one-dimensional model of Krticka & Kubat H2001fl . 

The DAC profile is calculated using the solution of the 
radiative transfer equation in the absence of any source 
function term in a point-star approximation as 



T{v,cj)) 



dr, 



(14) 

where n is line-opacity (its arbitrary value is selected to 
obtain optical depths of order unity) and vq descibes pro- 
file broadening as a sum of temperature broadening and 
microturbulent broadening by 



— ■'^th + ^turbJ 



(15) 



where Vturb is a microturbulent velocity. 

Surprisingly, the moving absorption features (DACs) 
in the spectra calculated from the hydrodynamical simula- 
tions are not predominantly caused by density variations. 
This can be seen from the fact that the wind optical depth 
can be approximated by the so-called Sobolev approxima- 
tion in the form of 

r~4. (16) 



Even with relatively small variations of wind density and 
velocity the variations of velocity gradient and also of the 
Sobolev optical depth may be large. This is especially the 
case when a typical feature of our hydrodynamical cal- 
culations develops, the so-called Abbott kink. In the re- 
gion of an Abbott kink the velocity derivative changes its 
sign from positive to negative and back. Thus, as can be 
seen from Eq. H16() , the optical depth variations due to the 
Abbott kink may be huge. 

3.2. Kinematic model inversion 

For the inversion of the DAC data (calculated using hydro- 
dynamical simulations) we use essentially the same code 
as Brown et al. H200H|I . The DAC data velocities are nor- 
malised with respect to the maximal velocity obtained in 
the hydrodynamical simulations, and only optical depth 
variation data for lower velocities are used in the kine- 
matic model inversion procedure. Finally, for the inver- 
sion of the velocity law, we take into account that the 
DAC profiles are calculated using data from only finite 
intervals of radius. 

4. Results 

For our simulated study we selected the parameters of 
the star HD 64760 because its wind variability was exten- 
sively studied in the course of the lUE "MEGA" campaign 
(Prinja et al. 11995(1 . Stellar and wind parameters used 
for the hydrodynamical simulations are given in Tabled 
Stellar parameters were taken from Hoffleit & Jaschek 
(HHHl, Humphreys & McElroy lfTUM|l . and Schaller et al. 
P992' see Fullerton et al. 11997(1 and for the parametriza- 
tion of radiative force we used generic parameters obtained 
by Pauldrach et al. The value of the microturbu- 

lent velocity was selected to obtain smooth wind profiles. 
Its value does not much influence the results obtained. 

Comparisons of data from hydrodynamical simulations 
and from inverted data are given in Figsm^jfor spots with 
relative amplitudes A = 0.05, 0.2, 0.5. We will discuss the 
reliability of each inverted relation separately. 

Table 1. Stellar and wind parameters of HD 64760 



Parameter 


Value 


R* 


22i?0 


Kot 


238 kms"^ 




23 100 K 


M 


20 Mq 


L 


1.210^ 1/0 


a 


0.650 


k 


0.050 


3 


0.28 


"turb 


50 kms"^ 



Apparently, the inversion of the mean velocity law 
Vr{r) is very good (middle left panels). Clearly, because 
during the inversion process we assumed that the radial 
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onset of Abbott kink (simulotion) 
end of Abbott kink (simulotion) 




F(r.,p)/F..(r) 



nversicn 

fluid parcel (simulation) 
Abbott kink (simulotion) 




0/2^ 

Fig. 1. Comparison of hydrodynamical simulations and inverted data for a spot with relative amplitude A = 0.05. 
Upper left panel: Variation of DAC optical depth (calculated as a difference to minimal absorption template) derived 
using the data from simulations. Upper right panel: Inverted variations of effective density and the locations of the 
onset and end of Abbott kink. Middle left panel: Radial velocity Vr{r, (f) taken from simulations at 18 equally spaced 
azimuths and the inverted mean radial velocity Vr{r). The Abbott kink in the hydrodynamical calculations can be 
clearly recognized. Middle right panel: Radial variation of azimuths of a fluid parcel relative to the spot, inverted 
location of density peak in azimuth and the azimuth of the onset of Abbott kink. Lower panels: Variation of density 
from the simulations p{r^ <j)) relative to the the density of ID flow without a spot /9id(^) and inverted density variation 
i^(r, (/)) for two different radii. 
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Fig. 2. The same as Fig^ however for spot with relative ampUtude A — 0.2. 



wind velocity is independent of azimuth, it cannot describe 
the azimuthal variations of wind velocity caused by the 
outward motion of an Abbott kink. However, it is possible 
to obtain a very good inference of the 0-averaged velocity 
law, at least for low-amplitude perturbations. The main 
reason is that the inversion method used is able to find the 



velocity law even for a nonvariable wind via the inversion 
of the unperturbed absorption profile. 

On the other hand the interpretation of the inverted 
rotational law is more complicated. The optical depth vari- 
ations of the line profiles calculated using hydrodynamical 
simulations are caused mainly by the velocity variations 
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Fig. 3. The same as Fig^ however for spot with relative ampUtude A = 0.5. 



due to the Abbott kink. Thus, we might expect that the 
inversion method would trace the movement of the Abbott 
kink and not the movement of real fluid parcels. To test 
this consideration, we calculated the radial variation of 
the azimuth of onset for the Abbott kink relative to the 
stellar spot. We define the onset of the Abbott kink as 



the point with the minimal radius where dv^jdr < 0. We 
have compared the radial variations of the azimuth of the 
onset for the Abbott kink and of the azimuth of the fluid 
parcel (calculated for constant — 0.4; results for other 
other azimuths are very similar) with the radial variation 
infered from inversion. These results are also plotted in 
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the upper-right panels of Figs. Inspecting Figs. ^El 
we conclude that the 'density' locus A(j){r) found by the 
inversion method in reality always traces the location of 
an Abbott kink and in fact does so with a high degree 
of reliability in those parts of the stellar wind where the 
Abbott kink is already formed. 

We plotted also (lower panels of FigsQllSl) the rela- 
tive azimuthal variation of the relative CIR density pro- 
file F{r, (f) obtained using the kinematic inversion method 
and the actual density profile calculated using hydrody- 
namical simulations (plotted relatively to the unperturbed 
wind model). Their comparison shows that the inversion 
method is able to infer density profiles well only relatively 
near to the stellar surface, where the optical depth vari- 
ations are not yet dominated by the variations of the ve- 
locity gradient. This is not the case in the outer parts of 
the wind where both density profiles are clearly different. 
This is caused by the formation of the Abbott kink, which 
enhances the optical depth. Note however that the spot az- 
imuth and its width at the stellar surface can be relatively 
well inferred from the kinematic inversion method. 

The fact that the main contribution to the DAC op- 
tical depth is given by the moving Abbott kink also does 
not allow us to invert to get correct wind densities. As dis- 
cussed in Sect.Elthe original method of Brown et al. (2003 ) 
assumes that the relative variations of wind density pro- 
file in azimuth are the same throughout the stellar wind. 
This also imposes limitations on the azimuthal variations 
of DAC optical depths. The situation is different with 
the azimuthal optical depth variations calculated with the 
data from simulations, since the main contribution to the 
DAC optical depth in this case is given by the moving 
Abbott kink. However, it is possible to slightly modify 
the original Brown et al. H2()()3|l method to obtain the ef- 
fective wind density that produces the DAC from simula- 
tions. Clearly, this effective density does not trace the real 
wind matter density, but rather the optical depth varia- 
tions, which are mainly given by the variation of velocity 
gradient due to the Abbott kink (see Ea. ljlGll 'l. This sit- 
uation is demonstrated in lower panels of Figs^OI where 
azimuthal density variations from simulations and from in- 
version are compared. Apparently, relatively near the stel- 
lar surface, the variations of velocity gradient are lower, 
and the optical depth variations arc influenced mostly by 
density variations. Thus, it is possible to infer relatively 
well the correct location of a footpoint of a spot at the 
stellar surface and its geometrical width. However, in the 
outer parts of the stellar wind, where the Abbott kink is 
already formed, the contribution of velocity gradient vari- 
ations to the optical depth variations is dominant. Thus, 
the inversion method in fact traces the movement of the 
Abbott kink. This can be seen from upper left panels of 
Figs ^"131 where density variations obtained from inversion 
are compared with the location of the onset and end of 
Abbott kink (defined as the radius where velocity gradi- 
ent changes its sign from positive to negative and from 
negative to positive respectively). 



5. Discussion 

We applied the kinematic inversion method of Brown et 
al. (P003) to DAC profiles calculated using hydrodynami- 
cal CIR models. We compared CIR structure inferred by 
kinematic inversion with the actual simulated one, as an 
important test of reliability of the inversion method. 

Although the kinematic inversion method is not able 
to calculate azimuthal variation of radial velocities, we 
showed that it gives reliable results for mean wind radial 
velocities even for relatively bright surface spots. However, 
because the highest optical depths of DACs in the nu- 
merical simulations arise mainly from spiralling velocity 
plateaux (Abbott kinks) and not from density enhance- 
ments, the inversion method does not trace the movement 
of real fiuid parcels, but of an Abbott kink. Thus, it is 
possible from DAC data to trace the spiral law of Abbott 
kinks propagating in the perturbed stellar wind. Near the 
stellar surface the spiral law of these patterns is inferred 
with a high degree of accuracy. This is with quantitative 
accuracy and gives a general indication of where absorbing 
matter lies. Finally, it is not possible to invert to get a reli- 
able density profile of the CIR, due to the above mentioned 
dependency of optical depth variations on the occurence of 
velocity plateaux, and due to the non-monotonic velocity 
law. Clearly, density variations are always overestimated. 
The density profile has the correct amplitude only very 
near to the stellar surface, where the velocity gradient 
variations do not significantly influence the optical depth. 
However, it is possible to obtain the correct azimuth of 
the CIR footpoint and its extent at the stellar surface. 

Note that we do not require that the wind variability 
manifested by the presence of DACs in the wind profiles 
of hot stars be caused by hot spots at the stellar sur- 
face. The hot spot model provides merely a convenient 
base perturbation for the generation of CIRs. The DAC 
profile calculated using this CIR model predicts all basic 
properties of observed DACs. Thus, it provides an appro- 
priate cornerstone for any technique aiming to infer real 
properties of the material causing DACs. We can easily 
test which properties of CIRs can be inferred from the ob- 
served DAC spectra and which cannot. We conclude that 
kinematic inversion method is able to infer some useful in- 
formation on several properties of the material causing the 
DACs and thus can be potentially used for observational 
tests of future DAC theories. 
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